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Dual representation of Polyakov loop 
in 3d SU(2) lattice Yang-Mills theory 

Florian ConradjB 
Institute for Gravitational Physics and Geometry, Physics Department, 
Penn State University, University Park, Pennsylvania, U.S. A 

We consider the expectation value of a Polyakov loop in 3d SU(2) lattice Yang-Mills theory and 
transform it to the dual representation in terms of sums over spins. The spin dependence of the 
amplitudes is computed explicitly by a graphical method. We also determine the asymptotic (large 
spin) limit of the amplitude factors. 



in 



^ 1 1. INTRO DUCT CON 

Strong-coupling expansions of lattice gauge theory provide an appealing physical model for the large-distance part 
of the quark potential [l[ : the sum over strong-coupling graphs corresponds to fluctuating strings of flux and naturally 
explains the confining potential and its iV-ality dependence. The problem is that these sums are hard to analyze in 
the continuum limit: they are either not convergent at weak coupling, or they are convergent, but difficult to evaluate. 
As the coupling decreases, a growing number of increasingly complex graphs has to be summed. 

Here, we will deal with the latter type of strong-coupling expansion, i.e. the one which converges for strong and 
weak coupling, ft results from an expansion of plaquette actions into a basis of characters, and from a subsequent 
integration over the connection 0, Thus, the sum over graphs is not an expansion in powers of /3, but rather a dual 
representation that is equivalent to the original lattice gauge theory [J, |5|, [6|, [7| . It can be viewed as a resummation of 
CV an expansion in /3. For this reason, we try to avoid the adjective "strong-coupling" and call the graphs instead spin 
foams 1). Originally, this name was introduced for SU(2) @, but it is also used for general gauge groups. In the case 
of SU(2), one obtains a sum over spin assignments to the lattice that satisfy certain spin coupling conditions. Each 
admissible configuration is a spin foam. 

In the early literature on the strong-coupling expansion, weighting factors were only determined for spin foams up 
to a certain complexity. Later, amplitudes were analyzed in general, and given as a function of arbitrary spin foams. 
For the partition function of 3d SU(2) lattice Yang-Mills theory, the form of the complete sum was first determined 
by Anishetty, Cheluvaraja, Sharatchandra and Mathur [J]. For the same theory, Diakonov and Petrov computed in 
• some detail the dual transform of a Wilson loop [6]. In an abstract form, Oeckl and Pfeiffer wrote down the dual 
representation for any dimension d>2, any compact gauge group G and for general observables [7|. 

In this paper, we present two new results in this direction: we use a graphical scheme to determine, in complete 
detail, the dual transform for a Polaykov loop in 3d SU(2) Yang-Mills theory. The graphical method makes the 
derivation more transparent and easier to check than a purely algebraic calculation. This result is used in ref. Q to 
JL" ' derive an exact string representation for two Polyakov loops. 
. £h ! The second result concerns the asymptotic limit of the amplitudes for large spins. The large spin limit is important 
^ for understanding weakly coupled processes in the dual representation. As was suggested by Diakonov & Petrov p, 
and further explored by the author [To| , dual gluons arise as spin waves in a weak-coupling approximation. This was 
shown by an argument that employs the large spin limit of 6j-symbols. 

In the presence of a Wilson loop, the argument is incomplete, however, since the amplitudes involve also 9j-symbols, 
and for these the asymptotic behaviour is not known so far [IfJ]. We improve this situation with the second result of 
the paper: we determine the large spin limit of the amplitudes for the Polyakov loop. This is possible, since we have 
chosen a zig-zag path for the loop, for which the amplitudes factorize into 6j-symbols. 

The paper is organized as follows: in section[TTl we specify 3d SU(2) lattice Yang-Mills theory with the heat kernel 
action. In section [LH1 and HVl we describe the dual transform of the partition function and of the expectation value of 
a Polyakov. The derivation is given in the appendix. 



43 



'Electronic address: fconrady@perimeterinstitute.ca 



Figure 1: (a) Triangulation T that is obtained by adding diagonals to cubes of k*. (b) Modified triangulation T with double 
edges where the Polyakov line passes through. 



II. SU(2) LATTICE YANG-MILLS THEORY IN 3 DIMENSIONS 

The partition function of 3-dimensional SU(2) lattice Yang-Mills theory is defined by a path integral over SU(2)- 
valued link (or edge) variables U e on a cubic lattice n with periodic boundary conditions: 



Z = J ( e n dt/ e ) exp ( - Sf(U f ) 



(1) 



The face (or plaquette) action Sf depends on the holonomy Iff around the face. We choose Sf to be the heat kernel 
action (for more details on the definition, see [HI). The heat kernel action has a particularly simple expansion in 
terms of characters, namely, 

exp (-SHE//)) =^(2 J + l)e-f^ +1 ) Xj (C/ / ). (2) 



The coupling factor (3 is related to the gauge coupling g via 

4 1 
ag z 3 

The expectation value of a Polyakov loop C in the representation j is 



P = — + o ■ ( 3 ) 



/ ( ndtk) ^J U c exp ( -J2 s f( u f)) ■ 



faU c ) = / ( EI dU e ) ttjUc exp ( - > j S f (U f ) ) . ! 4 ) 

Uc denotes the holonomy along the loop C. 



III. DUAL REPRESENTATION OF PARTITION FUNCTION 



Let us split the set of cubes into two subsets, called white and black, or even and odd, so that we obtain a 3- 
dimensional "checkerboard". Then, go to the dual lattice k* , and call vertices even if they are dual to even cubes, and 
odd if they are dual to odd cubes. Connect all odd vertices by edges. The resulting new complex is a triangulation 
T (see Fig.UJi). 

With these conventions, the dual transform of the partition function takes the following form: 

z = E + (uA ( II (-lyw-e-f-w-W] (5) 

{j c } T VeCT / VtCT / VeCK* / 

Each configuration {j e }r is an assignment of spins j e to edges e of T such that for each triangle of T the spins satisfy 
the triangle inequality. The edges of T belong to two groups: edges that are identified with edges in k*, and diagonal 
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Figure 2: (a) Zig-zag path of a Polyakov loop through the lattice ft. Only even cubes are shown, (b) The Polyakov loop in a 
2d slice of the lattice k. The arrows indicate how lattice points are identified. 

edges that were added to k* in order to form the triangulation T. In the amplitude, every edge contributes with 
the dimension 2j e + 1 of the representaiton j e . In addition, edges of k* give a sign factor and an exponential of the 
Casimir of the representation. For each tetrahedron, we get an amplitude factor 



3i 32 33 
34 35 36 



(6) 



where the spins jx, j 2 and j 3 are read off from any triangle in the tetrahedron. The spins j'4, 35, and jg are the spins 
on the edges opposing those of j% , 32 and j 3 . 

For large spins, 6j-symbols are approximated by the Ponzano-Regge formula [l2j 

A t «-i=(e^+e-^) . (7) 

Here, V t is the volume of the tetrahedron t when its edges have length j e + 1/2, and R t is its contribution to the 
Regge action: 

flt = X> + 1 /2)0 te + \ (8) 

eCt 



If we assume that large spins dominate, we can use this to obtain an asymptotic version of the spin foam sum: 

(9) 

We sum over signs s t = ±1 for each tetrahedron, due to the two terms of opposite phase in formula ijTj). 



z = e Efri(^+ 1 )) fn^^(^))««pf- i E 2 »-j--E|^.+i)' 

{je}r{st} VeCT / \tcT V 4 / \ eCfc« eCu* ^ J 



IV. DUAL REPRESENTATION OF POLYAKOV LOOP 

We now come to the dual transform of a Polyakov loop. We choose a Polyakov loop that runs along a zig-zig path 
through the lattice, and adopt suitable boundary conditions (see Fig. [2^ and b). More precisely, we assume that the 
Polyakov loop runs through a 1-2-plane of the lattice, and that the sites of the lattice are given by points 

(V2(k,l),m) 

(xi,x 2 ,x 3 ) = < or for (k,l,m) € Z 3 (10) 

t (x/2(fc-±,Z-i),m) 

with boundary conditions 

(xi + V2L,X2,x 3 ) ~ (xi,x 2 + V2L,x 3 ) ~ (x u x 2 , x 3 + L) ~ (xi,x 2 ,x 3 ). (11) 



Figure 3: (a) Extended tetrahedron that results from adding an edge of spin j to a degenerate tetrahedron of f. (b) Two 
consecutive degenerate tetrahedra t\ and t2 in the modififed triangulation T. 

That is, in each 1-2-plane, we identify points on opposing ends of diagonals, and in the cC3-direction we use ordinary 
periodic boundary conditions. 

The zig-zag path is not essential for determining the dual transform, but it helps us to find the asymptotic limit of 
the amplitudes. If we used straight lines, the amplitudes would contain 9j-amplitudes, and for these the asymptotic 
form is not known so far. For zig-zag paths, on the other hand, the amplitudes reduce to 6j-symbols and we can use 
the formulas by Ponzano & Regge and Edmonds. 

To describe the sum over spin foams, we have to modify the triangulation T: consider all faces of the dual lattice 
k* which are dual to edges of the Polyakov loop. In each of these faces, we add a second diagonal edge. The resulting 
complex is the modified triangulation T (see Fig. [lb). 

Then, the dual representation can be specified as follows: 

(t^c)=| e (n(^+ 1 )) (n^j (n (-^e-^^A . ^ 

{je}f \eCT / \tcf j \eCre* / 

Each configuration {j e }f is an assignment of spins e to edges of T such that 1. for each triangle of T the spins satisfy 
the triangle inequality, and 2. at every double edge the spins satisfy the inequality 

\ji-j\<j'i<ji+j. (13) 

Compared to the amplitude in ([5]), there is only one difference: in every cube of the dual lattice k* where the 
Polyakov loop passes through, we get an additional amplitude factor A t . To understand this, consider Fig.[Tj3: here, 
the Polyakov loop enters through the front of the cube, and exits on the right side, so we have added two edges — an 
edge with spin j' 3 on the front and an edge with spin j[ on the right side. Due to the presence of these two edges, we 
can define an additional degenerate tetrahedron in the cube. Apart from the tetrahedron defined by the edges ji, fa, 
j3, j4, j5 and jg, we get a tetrahedron formed by the edges j±, js, ]2, j'3 and j[, where the sixth edge is shrunk to 
zero length. 

Thus, we have two tetrahedra in the middle of the cube: the one we already had in Fig.[T^, and another, degenerate 
one. For this second tetrahedron, we receive an additional amplitude factor A t which is defined as follows: take the 
degenerate tetrahdron, and extend it by inserting an edge of spin j at the corner where the two double edges meet 
(see Fig. [3^). To this tetrahedron, we associate the factor 

A t = (-1)31-4 (_l)ii+ia+i3+j I h 33 32 1 _ ( W ) 

{ 33 Ji 3 J 

The spins j%, ji and j'3 are read off from one of the two triangles not containing j: if the short edge with spin j is 
drawn at the top (as in Fig. [3^), this triangle is on the left side of j in the direction of passage of the Polyakov loop, 
i.e. on the left side in the direction from j'3, j' 3 towards j\, j[. 

With this additional rule, the spin foam sum (fP2|) is completely specified. As in the case of the partition function, 
we can derive an asymptotic representation, assuming that all spins are large except the one of the Polyakov loop. 
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To deal with the asymptotics of the degenerate tetrahedra, we use Edmonds' formula jl3j : 

h h 32 X _ (-1 



f 3 3i 3 J A /(2j 1 + l)(2j 3 + l) 



dlM , (15) 



d J denotes the Wigner <i-function in the representation j and the angle i?t is defined as the angle between the two 
edges connected to j in the extended tetrahedron (see Fig.[3^i): 

q h(h + 1) + 33(33 + 1) - h(h + 1) f , - 

cosv t = = (16) 

2 v / ii(ii + i)i3(i 3 + i) 

Let us apply this formula to two consecutive degenerate tetrahedra in T (see Fig. Et"). They give the amplitudes 

/_iy'i+j' 2 +fa+j J 3i h h 1 f]\ji+hi+hb+j J h 3u hp 1 ^) 
i 32 ji 3 J i 3u 3i 3 J 

1 j 1 ' \h + k h + l j j [ L > \hi + m h + l 3 

= { - iri 7mTm^ diM TmMmm ^ (18) 

= v^ 2 + Wi T dv^u + 1 ( - 1)2fc+m " ( - ir "' ^ (19) 

' 4^4n(^)- (20) 



V2j 2 + 1(2 ji + l)V2ju + 1 



If we multiply these factors for the complete sequence of degenerate tetrahedra U along the Polyakov line, the total 
result becomes 

<tr^c>*4 E zfn^+i)) (\i i:/ H +] .,■ ( 2i ) 
n ^7^w exp ( istiit ) U^f-iE 2 ^- eJ^^ +i )) ■ ( 22 ) 

^non — degenerate tcf / ^ eCft* eCu* / 

In this expression, the sequence of double edges is numbered by i — 1, . . . , N. The associated spins are denoted by ji 
and 2% and "&i is the angle between edges number i and i + 1. 



V. DISCUSSION 



In this paper, we computed the explicit dual transform of the expectation value of a Polyakov loop in 3d SU(2) 
lattice Yang-Mills theory. To improve the transparency of the derivation, most calculations were done in a graphical 
scheme (see appendix). In ref. [9|, our result is used to derive an exact string representation of two Polyakov loops. 
Due to the specific choice of the loop — a zig-zag path — the amplitudes factorized into 6j-symbols and we could 
determine their asymptotic large spin limit. Such a limit is useful when one tries to understand how spin waves and 
Coulomb force emerge in the dual representation (see ref. 10]). Our asymptotic expression may be also helpful for 
analyzing the sign behaviour of the amplitude, which is an important issue in Monte Carlo simulations [bH ]. 

As a side remark, we point out a similarity between the asymptotic formula (|22p and the amplitude for spinning 
particles in 3d quantum gravity (see eq. (153) in appendix B.2 of [ill). In formula ([22)l . the Polyakov loop appears as a 
product of Wigner d-symbols in the representation j . A similar structure is obtained when we consider a loop-shaped 
Feyman diagram of a spinning particle in 3d quantum gravity. Such a relation is not entirely surprising, since 3d 
SU(2) Yang-Mills theory can be regarded as a deformation of 3d quantum gravity (see e.g. jldl. H^V 
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Appendix A: GRAPHICAL DERIVATION OF DUAL REPRESENTATION 

In this appendix we derive the dual representation of the expectation value ([J) by a graphical method. First we 
will explain our graphical notation and express a number of identities in terms of it. Then, we apply these identities 
to obtain the dual representation. 

1. Conventions 

There exist various graphical schemes for representating calculations with SU(2) tensors. A standard refererence 
is the work by El-Baz & Castel [17} : it incorporates the earlier notation by Yutsis, Levinson & Vanagas and Brink 
& Satchler pilOjj], and extends it by symbols for representation matrices and group integrals. Our present choice of 
notation will be similar to these references as far as invariant tensors are concerned. For representation matrices and 
integrals, we prefer to use a different convention which is inspired by ref. [20l |. 

A directed line with label j represents the identity in the spin j representation of SU(2): 



3 



(Al) 



Similarly, a representation matrix is symbolized by 



U 

— •- 



3 
-4- 



D j (uy 



(A2) 



The target of the line corresponds to a vector in the j-representation Vj , while its source corresponds to a one-form 
in the dual V* . The dot indicates the dependence on the group element U, and the arrow on it specifies whether we 
have Di(U) or Z? J (C/ _1 ). Integrations are indicated by "cables" around one or several lines: 



m 





h 




— 4 — 


4 — 


— 4 


h 


4 — 










4 — 


— 4 


33 





mi 

m 3 



dU D^(U) m i mi D^(U) m - m2 D^{U) m ^ 



(A3) 



We use 3jm-symbols to construct normalized invariant tensors for triple tensor products of irreps. For such tensor 
products, the subspace of invariant tensors is one-dimensional, so we only need to fix one normalized tensor for each 
type of space. For the tensor product Vj 1 ® Vj 2 ® Vj 2 we take the 3jm-symbol itself, and symbolize it by a node with 
three outgoing arrows: 



mi 



3 1 



32 



m 2 



.73 



Jl 32 33 

mi m2 m3 



(A4) 



7713 



The order of spins is specified by a circle on the node. This is necessary, since permutations change the sign: 



Jl 32 33 

m\ mi ni3 



(-1) 



31+32+ J3 



31 33 32 

mi m3 vri2 



(A5) 
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For the dual space V* ®V* 2 ®V* 2 ,we choose the dual of [|S4)) : 
mi / mi 



3 1 



m 2 



32 



33 

-4 m 3 



3i 



\ 



m-2 



32 



33 

■> m 3 ) 



3i 32 33 
mi m 2 m3 



(A6) 



The components are again those of the 3j'm-symbol, since the 3jm— symbol is real. To construct normalized invariant 
tensors with upper and lower indices, we define the tensor 



-• 4- 



y/Wi + 1 x 



m 












t 




3 




f 







—4 



in 



(A7) 



It is a normalized and invariant tensor in V* ® Vj*. By contracting it with the tensor (|A4|) we obtain a normalized 
invariant tensor in Vj 1 (8 Vj 2 ® V^* : 



mi 



mi 



3i 



m 2 



.72 



J3 

-4 m 3 



m 2 



J2 



« 33 

-> • «- 



m 3 



E(-!) i£ 



3+"»3 



J m3,-mg 



3l 32 33 \ _ f]\fo+m 3 ( 3l 

mi m 2 m 3 J \ mi 



3i 32 33 

m 2 -m 3 



The invariant tensor in V* <g> V* (£> V,-, results from dualization. 

Jl J2 J J 



mi 



3i 




m 2 



J3 

* m 3 



/ "ii 
Ji 



V 




m 2 



J3 

^ m 3 / 



' mi m 2 — m 3 



The remaining possibilities are fixed by applying the same procedure to j± or j 2 instead of j 3 . 



2. Identities 



(A8) 



(A9) 



Next we express a number of identities in our graphical notation: we have that 

i 



U J A 
-• 



3a 



■ lb 

■ 3b 



Jb 



£/ 31 V 



(A10) 
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Figure 4: The loops from the character expansion are attributed to even cubes. Application of the identity (IA10[) leads to the 
diagram on the right-hand side. The loops are drawn without the dots from representation matrices. 



This is equivalent to eq. (13), p. 296, in Varshalovich et al. [21|. We also need the identity 



32 



33 




(A18) 



EM) 



ie+P+Ji+n+h+o 



3i 
k 



35 




36 

-p 



32 

I 



36 

p 



—n 



33 34 35 

m n —o 



3i 32 33 I / 3i 32 33 
34 is 3a ) \ k I m 



(A19) 



The latter is equivalent to eq. (6), p. 454, in [2l| . 



3. Derivation of dual representation 



The transformation to the dual representation proceeds in two steps. First we expand the plaquette action for each 
face into characters: 



exp(-S f (U f )) =£(2j> + 1) 



XiA u f) 



(A20) 



3f 



In the second step, we integrate over the connection variable U e , so that there remains only a sum over spin variables. 
Before and after the integration, we use, in addition, identity (jAlOj) and (1A18|) . This allows us to reduce all tensor 
contractions to 6j-symbols. 

As explained in the main part, we separate the cubes of the lattice k into alternating even and odd cubes. Each 
face of k is attributed to an even cube. In our graphical representation, the trace Xjf(Uf) in (|A20j) becomes a loop 
of spin jf. Thus, we have a loop of spin jf for each face / of an even cube (see Fig. [4]). On each pair of neighbouring 
strands we can apply identity (|A10[) and thereby turn the six loops into a single diagram. For each edge of the cube, 
we get a new strand, a sum over its spin, and the dimension of the representation. 

To proceed further, we need to label the spin variables. We do this by labelling the vertices, edges and faces of the 
even cube and its immediate vicinity (see Fig. O . As in Diakonov & Petrov's paper 0] , faces are labelled by capital 
letter A, B,C . . ., edges are given numbers 1,2,3..., and vertices receive lowercase letters a,b,c. .. The spins of loops 
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Figure 5: Labelling of faces, edges and vertices of even cubes in «*. 

are denoted by Jx where X is the label of the face to which the loop belongs. For each edge i, we get two new spins, 
coming from two even cubes that share the edge: we call them ji and j[. Magnetic quantum numbers are associated 
to vertices and designated by letters o,p,q . . . with a subscript for the relevant vertex. 
Thus, we obtain the following, labelled diagram for the even cube: 



n ^ + x ) 



X 12 



o 




(A21) 



Strands along edges carry spins of type ji, while the strands in the corners have the spin Jx of the corresponding 
face. Identity (|A10[) requires that nodes on opposite ends have opposite orientation. We satisfy this by choosing a 
common orientation for all three nodes near a vertex and by letting this orientation alternate as we go from corner 
to corner. In this manner, we receive two kinds of vertices: vertices of type I that have clockwise orientation from 
outside the cube (a, c, f,h), and vertices of type II with anti-clockwise orientation from outside the cube (&, d, g, e). 
We indicate these orientations by eight circles: four of them are situated in front of the cube, and four of them lie 
behind it. This pattern of orientations is consistently extended to all even cubes. Then, a given vertex of k is either 
of type I for all surrounding even cubes, or of type II. 

Next we apply identity (|A18[) to replace the corners of the diagram (|A21|) by 6j-symbols and 3jm-symbols. For a 
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type I vertex such as a, we get 




(A22) 



(A23) 



For a type II vertex like b, we have, on the other hand, 




We repeat this procedure everywhere in the lattice and obtain a 6j— and 3jm— for every corner of an even cube. 
Altogether this gives us four 6j— and 3jrn-symbols for each vertex. For each edge, we receive, in addition, a sign 
factor and a factor (2jj + l)(2f z + 1). 

The 3jm-symbols at a vertex are connected to those at other vertices by representation matrices (see Fig. [6]). We 
integrate over these representation matrices using identity (1A12|I and (|Al4j) . The result is a contraction of tensors for 
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Figure 6: Contraction of 3j'm-symbols: the dots indicate representation matrices Dj(Ui). If the Polaykov loop passes through 
the vertices a and b, there is an additional strand of spin j (drawn in the middle between 3j'm-symbols) . 

each vertex. When a and b are two vertices in the Polyakov loop C, we obtain 




14 14' 17 



I II 

For every edge outside C, the integration forces fa and j[ to be the same and gives, in addition, a factor (2ji + 1) _1 . 

It is sufficient to evaluate these two spin networks, since the path of the Polyakov loop is periodic. The spin 
networks outside the Polyakov loop are the same with the edges of spin j deleted. As we see below, the spin network 
at a (type I) is related to that at b (type II) by a number of simple operations. Thus, we only need to compute the 
spin network of type II. We use identity (|A18[) to split it into two diagrams with 6 spins 1 : 




(A26) 



H 2 Hi 



1 This simple splitting of diagrams is possible, since we have chosen a zig-zag path for the Polyakov loop. For a straight Wilson line, we 
would receive a 9j— symbol. In that case, the reduction to 6j'-symbols is more complicated, and the asymptotic behaviour not known. 
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The modulus of these diagrams is equal to 6j-symbols. It remains to compute their phase. We do this by translating 
the graphs back into formulas. For that purpose, we have labelled each edge by lowercase letters to indicate the 
magnetic quantum number of the 3jm-symbols. For diagram Hi, we receive 

Hi = (-l) j9+Vb 5 v , - w A-iy o+w " ( j>1 j ' 2 jg \ (-l) jl7+r "( i' 1 jl7 Jl8 
v ; Vb, w b \ j \u b v b -w b J v ' \u b -r b x b 

x (-iY i »+ x b ( i' 2 Jl8 • ?16 ) (_iy'ie+«& i 3s he 3n \ 
\v b -x b s b J \y b -s b r b J 

_ (_ifjo ( fi 32 h \ /_iyi8+x b ( fi 317 318 
\ -u b -Vb -yb J V ~ Ub ~ Vb Xb 

X (-iy™+ s » ( J ' 2 JlS Jl6 \ (-l)n7+r b ( 39 316 317 

\-v b -Xb s b J y -y b -Sb n 

2j« ( 3l 32 39 \ (_-i\ha+x b ( 3l 317 318 



1 u b v b y b J \u b r b -Xb 



\ji6+s„ ( 32 318 316 \ /_,y 17 +rs ( 39 jl6 jl7 



V v b x b -s b J \Vb s b -r b 

C_]_)2j 9 +E 4 2ji / 3l 32 39 1 _ /_ 1 x2j 9 +2^+2ji 7 j j'l 32 39 
\ 316 317 318 J I 316 317 318 



C_^2ji+2ji 7 J 3i 32 39 I (A27) 

' 316 317 318 ' 



Diagram II 2 gives 



II _ r_iy'i+39+j' 2 ( 3i 32 39 \ i^y+mt ( 3i j'i 3 \ 
' \ o b p b w b J v ' \o b u b -m b J 



32 3 32 \ (_->y s +w b I 39 32 3i 



x JZ J J2 (-1) 

\p b m b v b J y \ -w b v b u b 

f_iy'i+h+j' 2 ( 3i k 39 \ t]_y+m b ( 3i 3\ 3 
\ o b p b w b J \ o b u b -n 



x [ 32 j j'2 \ r_iy 9+v>b ( 3 9 32 fi 
\ p b m b -Vb J V -Wb -Vb u b 

= (— \y'2+ v t>+iv+ w b+fi+ u b ( 3i 32 39 \ r_^\j+m b ( 3i 3\ 3 

\o b p b w b J v ' \o b u b -n 

x C_iy'2+«6 ( 32 j 32 \ iiy'x+ub ( 39 32 3i \ 
V ' \Pb m b -Vb J \ -w b -Vb u b J 

= (_iy'2+j9+j[+2u b | 3i 32 J9 1 _ (_iy' 1 +j' 2 +j 9 +2j' 1 I j'i j'2 39 1 (A28) 
V ' 1 j'2 j'l j J 1 32 Jl J J K ' 
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Together this yields 



'A 


32 


39 \ 
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1 J 






39 


1 / fl 32 39 


\32 


3i 




J \ jl6 317 318 



II = IIoIIi = r_iyi+Ji+JB+ a Ji <^ J1 J2 Jy K (-1) 2 H+Vi7 J Jl ^2 ja l / A29 ) 

' ' ' ' ' ' iie jn 3i8 J 

(A30) 

The spin network of vertex a (type I) is related to that of vertex b (type II) by the following steps: 

1. rotate the lower right corner of spin network II onto the upper left corner and identify its edges with those of 
spin network I, 

2. take the complex conjugate (i.e. reverse all arrows on edges), 

3. reverse the arrow on the node of edge 15, 

4. reverse the arrow on the edge j. 

Step 1. corresponds to a relabelling of spins in (|A30j) . Step 2. does not change the value, since we deal with real 
quantities. In step 3. we get a factor (— l) 2 -? 15 , since 

its ji B \ = ( _ 1)2jl5 ( jib j 15 \ (A31) 
w a y a J v ' \ y a w a J v ' 

while step 4. does not produce any factors: 

E(_]\j+ma | 3ii 3 Jl4 ] c_^y'+ma f Jl J 3i 
\ Pa -m a W a / Ua ~™a O a 



/ J14 3 Jl4 Jl J Jl 

^ 1 Pa ™a v a J \ u a m a o a 



Therefore, 





3' 2 


39 \ 
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39 


{32 


3i 


3 J 


{ 316 


Jl7 


318 



For vertices outside the Polyakov loop, this reduces to 



(A32) 



I = (-l) 2jI5 Hrelabelled (A33) 

= (-l) 2j15 f_iy'i4+Ji+il5+2ji3 / Jl J15 1 f J14 Jl J15 1 , A34 n 

\ fl 3U J J I J4 J13 J12 J 

= (-l)Ji+Ji4+Ji5+ 2 J4 / J' 1 315 1 f jl J14 J15 j / A35 n 

1 Jl4 fl J J I Jl3 J4 Jl2 J 

In summary: along the Polyakov loop, vertices of type I and II contribute 

J _ (_l)ji+jiA+h5+2ji | Ji Jl4 jl5 1 f jl jl4 jl5 1 (A36) 
1 Jl4 fl j J I J'l3 J4 J12 J ' 



(A37) 



IU= = 1 - =(-l) 2j 4 Jl J ' 14 J ' 15 }■ , (A38) 

U V2J14 + TVWl + 1 \J13 J4 J12 ' 

H| = 1 = /^r ( Jl J' 2 J9 I (A39) 

|J "° V2jl + V2j2 + 1 \J16 Jl7 J18 ' 
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Since vertices of type I and II alternate along the lattice, the factors of (— l) 2j4 and (— l) 2j17 cancel each other, when 
we multiply all vertex factors. Outside the Polyakov loop, the \/2ji + 1 factors combine to give (2ji + 1) . 

Let us summarize what we have obtained so far: the transition of Fig. |4] introduced new spins ji and j- and factors 
(2ji + l)(2j z ' + 1) for each edge. From the application of identity (|A18[) to diagram (|A21[) , we got a factor 



3k Jl Jrn 

Jx Jy Jz 



(A40) 



for each even cube near a vertex. Due to type II vertices, we also had a sign factors (— 1)^+^ per edge. The integration 
over group variables produced factors 





3l4 


315 1 




314 


315 1 






3 J 


{ 313 


34 


312 J 



I — (_iyi+m+n5 ) Ji- J\a ( ) jlo ( (A41) 



XI = (_i)Ji+J2+J9 ) Ji i i ji J2 jy i (A42) 

\ h 3i 3 J \ he 317 318 J ' 

for vertices along the Polyakov loop C, and outside of it 

I = / jl 314 Jl5 } , (A43) 

1 313 34 312 j 

II = ( 32 h ) . (A44) 

1 316 317 318 J 

For every edge outside C, the integration imposes ji = j[ and gives a factor (2ji + l) -1 . 

We now switch to the description in terms of the modified triangulation T: spins Jx on faces / C k turn into spins 
j e on edges e C n* , and spins ji and jl on edges become spins j e on diagonals in T. The spin foam sum takes the 
form 



(^Uc) = \ £ (n(2ie + i)) 

UJf \ecf / 



f n (-Hf n c-D^fn^fn 

\sijigle diag. e C f / \doublc diag. e, e' C T / ytcT / VCt* 



-|je(ie+l) f ( A45 ) 



where for tetrahedra as in Fig. [T^, 



and for degenerate tetrahedra (as in Fig. [TJ3) 



A, = < j; J » } , (A46) 



A, = j 3) J» « j. . (A47) 



Except for sign factors these formulas are identical to eq. (fT2"|) , ([5]) and (Ti"4"|) . 

To show equivalence with (TT2]), we need to shift the sign factors on diagonals to edges of n* . For this purpose, we 
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think of the triangulation as being made up of octahedra around vertices of k* : 




Let us first regard the case, where the Polyakov loop is absent. Then, T is the same as T and we have a sign 
factor (— l) 2je for each diagonal e of T. We can distribute these factors such that each factor belongs to exactly one 
octehadron. One possibility is, for example, that every octehadron carries the factors 



_1^2j 9 +2j 6 +2j 8 +2ji 2 +2j 7 +2j 3 _ (_-^2J A +2J F 



(-1) 



(A48) 



(-1) 



2J a +2J b +2Jc+2Jo+2Je+2J f 



(-1) 



2ji2+2jio 



(A49) 



When we multiply the octehadra, the factors (— l) 2 ji2+2ji cance i each other due to periodicity. Therefore, we have 
the identity 



2j e 



(A50) 



i.e. the sign factors can be shifted from diagonals of T to edges of n*. 

In the presence of the Polyakov loop, the argument is slightly modified. Depending on its location, the Polyakov 
loop affects edges of type 8 and 7 or 1 and 2. Suppose the Polyakov loop goes through 8 and 7. Then, we have factors 



(A51) 



from the double edges 8 and 7. We repeat our previous argument with the first factor on the right-hand side, and 
absorb the second factor into the tetrahedral amplitude. When the Polyakov loop passes through 1 and 2, nothing 
changes, since the factors of 1 and 2 are attributed to adjacent octahedrons. This gives us precisely the amplitude we 
described in the main part of the paper. 
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